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We discuss Standard Model (SM) and New Physics (NP) descriptions of D" mixing. The SM part 
of the discussion addresses both quark-level and hadron-level contributions. The NP part describes 
our recent works on the rate difference AFd and the mass difference AMd. In particular, we 
describe how the recent experimental determination of AMd is found to place tightened restrictions 
on parameter spaces for 17 of 21 NP models considered in a recent paper by Hewett, Pakvasa, Petrov 
and myself. 



I. INTRODUCTION 

Given the forthcoming operation of the LHC, per- 
haps the dominant role of experimental flavor studies 
in particle physics will be supplanted by discoveries in 
the so-called new physics. Even if flavor physics faces 
an unsure future, all would acknowledge its remark- 
able recent progress via the observation of rare phe- 
nomena such as CP-violation in i?-mesons or particle- 
antiparticle mixing for and mesons. If new 
physics is indeed observed, the continued exploration 
of rare observables could well be an asset in decipher- 
ing exotic LHC events. 

My purpose in this talk is to describe two recent 
theoretical contributions to mixing [l], Q . [1^ Both 
Ref. [3] and Ref. 2] should be considered in the con- 
text of the recent HFAG values Q, 
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In light of the Physical Review Letters criteria of 'ob- 
servation' (> 5a) or 'evidence' (Scr-to-Sa), we see that 
the above 2Aa determination for xu amounts to a 
'measurement' (< Scr). As such, we all await improv- 
ments in sensitivity for charm mixing. 

The observed signal is seen to occur at about the 1% 
level. Whether or not this is the magnitude expected 
for the SM signal is a topic I will discuss shortly. At 
any rate, I wish to also consider the possibility of a 
NP component in mixing amplitude, 

M,nix = MsM + A^NP . (2) 

The relative phase between A^sm and Mnp is not 
known. Thus, in our detailed study of various NP 
contributions to xd in Ref. 0] we most often compared 
the NP predictions to zblcr, ±2cr windows relative to 
the central xb value of Eq. H]) . 



/Lt = M to one at, say, the charm scale fi = rric. This 
takes the form 

{f\nNp\i)^Gj2G^ifi) {f\Q^\^)i^i) , (3) 



where the prefactor G has the dimension of inverse- 
squared mass, the Ci are dimensionless Wilson coeffi- 
cients, and the Qi are the effective operators. At the 
leading order of dimension six, it turns out that there 
are eight four-quark operators, 

Qi = (ul^^cl) {ulI^cl) , 

Q3 = (uLCfl.) (JirCl) , 

Qi = (u^cl) (uj?cl) , 

Q5 = {uro-^l^cl) {urct^^^cl) , 

Qe = (ur^j^cr) (ur^j'^cr) , 

Q7 = {ulcr) {ulcr) , 

Qs = iuLCTf^iyCR,) (uLcr^^CR) . (4) 

Any given NP contribution will often involve several 
of these, but in all events never more than these eight. 
The evolution is determined by solving the RG equa- 
tions obeyed by the Wilson coefficients. 



dlog/i 
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where 7 is the 8x8 anomalous dimension matrix [J]. 
The output of this calculation is a set of RG fac- 
tors ri(fj,,M) which are expressed in terms of ratios 
of QCD fine structure constants evaluated at different 
scales, e.g. as with 
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A. Operator Product Expansion (OPE) and 
Renormalization Group 

An important technical aspect of Refs. [H, 0] is the 
process of relating an amplitude at some NP scale 



B. Operator Matrix Elements 

One needs ultimately to evaluate the D^-to-D'^ ma- 
trix elements of the eight operators {Qi}. In general, 
eight non-perturbative parameters would need to be 
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evaluated by some means such as a lattice determi- 
nation. As a practical matter, the method used in 
Refs. [l|, [3| is to introduce a 'modified vacuum sat- 
uration' (MVS), where all such matrix elements are 
written in terms of the known matrix elements of (V- 
A)x (V-A) and (S-P)x (S+P) matrix elements Bu and 
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FIG. 1: Loop diagram for D° ^ D° 
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where the number of colors is Nc 
Ref. 0, we define 



B 



(S) 



B 



(S) 



(7) 

3 and, as in 
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With the above theoretical machinary in hand, we are 
now ready to consider SM and NP contributions to D'^ 
mixing. 



sin 9c, itself seems to suggest that the experimental 
signal (near the 0.01 level) is easily understood. But 
not so fast! For convenience, let us set rrid ~ 0. Then 
the only mass ratio that appears in the problem is 



0.006 



(9) 



Table U examines one of the loop-functions for APd 
and shows the results of carrying out an expansion 
in powers of z. We see that the contributions of the 
individual intermediate states in the mixing diagram 
are not intrinsically small - in fact, they begin to con- 
tribute at 0{z^). However, flavor cancellations re- 
move all contributions through ©(z^) for AFd, so the 
net result is 0{z^). Charm mixing clearly experiences 
a remarkable GIM suppression! 

We understand the reason for this. mixing van- 
ishes in the limit of exact SU(3) flavor symmtery. It 
is nonzero only because flavor SU(3) is broken, and 
indeed, D° mixing occurs at second order in SU(3) 
breaking [11 . A factor of z will accompany each or- 
der of SU(3) breaking and the rate difference t/d will 
experience an additional factor of z due to helicity 
suppression. 



II. STANDARD MODEL ANALYSIS 

One can use quarks or hadrons as the basic degrees 
of freedom in carrying out the SM analysis of mix- 
ing. In principle, these should give the same result. 
However, as we shall see, rather different features ap- 
pear in each description. 



TABLE L Flavor cancellations in AFd- 
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A. Quark-level Analysis 

At leading order in the SM, the OPE for D° mix- 
ing consists of two dimension-six four-quark opera- 
tors 01 ■ The next order contains fifteen dimension- 
nine six-quark operators. For each increasing order in 
the OPE, there are still more local quark and gluon 
operators and the problem of determining operator 
matrix elements becomes ever more severe. For this 
reason, the dimension six sector has received by far 
the most attention. 

The dimension six amplitude is depicted in Fig. [TJ 
Since the 6-quark is essentially decoupled due to the 
tiny Vub value, only the light d, s quarks propagate in 
the loop. The Cabibbo dependence of this diagram. 



Of course, this is just the leading order (LO) result 
in QCD, and we should consider the next-to-leading 
order result as well. 



_ (LO) (NLO) 
„ - „(LO) , „(NLO) 



(10) 



This has been done in Ref. [6| and the results are 
summarized in Table [Hi which reveals that ?/d is 
given by y^LO to a reasonable approximation (due 
to the removal of helicity suppression by virtual glu- 
ons) whereas xb is greatly affected by destructive in- 
terference between XhO and xnlo- The net effect is 
to render yj^ and xd of similar small magnitudes, at 
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TABLE II: Results at dimension-six in the OPE. 
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least through this order of analysis, as compared to 
the experiment signal. 

It is not inconceivable that the quark-level predic- 
tion of CCD and ?/d just described might be considerably 
affected by a higher order in the OPE Q which suffers 
less z suppression. Simple dimensional analysis 
suggests the magnitudes xd ^ yu ^ 10~^ might be 
achievable, although order-of-magnitude cancellations 
or enhancements are possible. 



are kinematically allowed. However, this cannot be 
true for four-pseudoscalars because decay into four- 
kaon states is kinematically forbidden. In Ref. [1] it 
is estimated that these 'kinematically-challenged' sec- 
tors can provide enough SU(3) violation to induce 
Ud ^ 10^^. I personally find such an argument to 
be an important advance in our understanding of the 
subject. At the same time, it is unfortunately more 
persuasive than compelling due to the uncontrollable 
uncertainties inherent in this line of reasoning. 

To summarize, we have just described how the ob- 
served D'-^ mixing signal could well arise from SM 
physics, but the associated numerical prediction is 
seen to be lacking in precision. This conceivably leaves 
room for some NP mechanism to co-contribute or even 
dominate the SM signal. In the following we consider 
in turn NP analyses of the width difference j/d and 
the mass difference xd- 



B. Hadron-level Analysis 

Most of the work involving the hadron dcgrce-of- 
freedom has been done on j/q. One starts with the 
following general expression for AFd, 



AFd = T^Ini/ 



(11) 



To utilize this relation, one inserts intermediate states 
between the |AC| = 1 weak hamiltonian densities 

n\^'^^=\ Although this can be done using either 
quark or hadron degrees of freedom, let us consider 
the latter here. Clearly, some knowledge of the ma- 
trix elements (n|7i|^'^'~^|Z?*') is required. 

One approach is to model | AC| — 1 decays theoret- 
ically and fit the various model parameters to charm 
decay data. Some time ago, AFd was determined in 
this manner and the result yn — 10~^ was found [ll| . 
This value is smaller than the recent BaBar and Belle 
central values. 

Alternatively, one can arrange for charm decay 
data to play a somewhat different role. The earli- 
est work in this regard focussed on the P+P^ = 
■K+n-,K+K-,K-TT+,K+Tr- states [H [ll. In the 
flavor SU(3) limit, this subset of states gives zero con- 
tribution due to cancellations. But SU(3) breaking 
had already been known to be significant in individ- 
ual charm decays. Since the study of charm decays in 
the 1980's lacked an abundance of data, these refer- 
ences could only conclude that '?/d might be large'. 

A modern version of this approach now exists, al- 
though the analysis takes an unexpected direction Q . 
Since SU(3) breaking occurs at second order in D° 
mixing, let us hypothesize that the contribution of the 
P+P" sector is in fact negligible due to flavor cancel- 
lations. Likewise for all other sectors whose decays 



III. NP AND THE WIDTH DIFFERENCE 

At first glance, it would appear unlikely that NP 
could affect yjj because the particles contributing to 
the loop amplitude of Fig. [T] must be on-shell. Since 
NP particles will be heavier than the charm mass, 
'there can be no NP contribution to j/d'- Or so goes 
the argument. 

However, as explained in Ref. NP effects in 

Tiilf^'"^ can generally contribute to yu- In the loop 
amplitude of Fig. [Jl the NP contribution (empirically 
small for AC = — 1 processes) arises from either of 
the two vertices. We represent the NP AC = — 1 
hamiltonian as (indices i,j,k,£ represent color), 

9,9' 
/ — 



Ci = UiTiqj qjT2Ci , 
©2 = UiTiql qjT2Cj , 



(12) 



where Dqqr and the spin matrices Fi_2 encode the NP 
model. Ci,2(Ai) are Wilson coefficients evaluated at 
energy scale fi and the flavor sums on q, q' extend over 
the d, s quarks. 

This leads to a prediction for the NP contribution 
to yi). For a generic NP interaction, one finds (with 
the number of colors Nc — 3) 

4:V2Gf 



2/D = 



MdFd 



5 

+ K2SuSjk)y] Ic.{x,x') {^\ Of'' (13) 



a=l 



where {Ka] are combinations of Wilson coefhcients, 
Ki = {CiCiNc + (C1C2 + C1C2)) , 
K2^C2C2 , (14) 
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and the {OJ^^} are four-quark operators written 
down in Ref . [T] . Numerical results for some NP mod- 
els are displayed in Table Hill 



TABLE III: Some NP Models and yu. 



Model 


Yd 


Comment 


RPV-SUSY 


6 ■ 10"^ 


Squark Exchange 




-4 • 10^'' 


Slepton Exchange 


Left-right 


-5 ■ 10^'' 


'Manifest' 




-8.8 ■ 10-^ 


'Nonmanifest' 


Multi-Higgs 


2 ■ 10-"' 


Charged Higgs 


Extra Quarks 


10^* 


Not Little Higgs 



One sees that the entries, aside from R-parity vi- 
olating SUSY, produce small contributions. We em- 
phasize, however, that Eq. p2|) and Eq. (|13p represent 
general formulae for the contribution of all NP mod- 
els of |AC| — 1 interactions, encompassing those not 
included in Table [ml 



IV. NP AND THE MASS DIFFERENCE 

As the operation of the LHC looms near, the num- 
ber of potentially viable NP models has never been 
greater. In this section, I will give an overview of 
Ref. 0, whose hallmark is the study of many (21 in 
all) NP models. Perhaps the best way to start is to 
consider the different ways that 'extras' can be added 
to the SM: 

• Extra gauge bosons (LR models, etc) 

• Extra scalars (multi-Higgs models, etc) 

• Extra fermions (little Higgs models, etc) 

• Extra dimensions (universal extra dims., etc) 

• Extra global symmetries (SUSY, etc) 

Although this approach does not provide a totally 
clean partition of NP models {e.g. obviously SUSY 
contains extra particles appearing in other categories), 
it proved useful to the authors of Ref. . 

The broad menu of NP models which were analyzed 
is listed in Table HVl The extensive content of this list 
{e.g. there are four different SUSY realizations and 
three involving large extra dimension) indicates how 
rich the field of NP models has become. Of course, 
the subject of NP is by now fairly mature (in prepar- 
ing this talk, I realized that my first paper on charged 
Higgs bosons |l4| was written nearly 30 years ago!) 
and thus many models have been well exposed to the 
scrutiny of experiment. This would seem to imply 
that parameter spaces for the various models have 
shrunk so much that a measurement like D'^ mixing 



TABLE IV: NP models studied in Ref. 0] 

Model 

Fourth Generation 
Q — —1/3 Singlet Quark 
Q = +2/3 Singlet Quark 
Little Higgs 
Generic Z' 
Family Symmetries 
Left-Right Symmetric 
Alternate Left-Right Symmetric 
Vector Leptoquark Bosons 
Flavor Conserving Two-Higgs-Doublet 

Flavor Changing Neutral Higgs 
FC Neutral Higgs (Cheng-Sher ansatz) 
Scalar Leptoquark Bosons 

Higgsless 
Universal Extra Dimensions 
Split Fermion 
Warped Geometries 
Minimal Supersymmetric Standard 
Supersymmetric Alignment 
Supersymmetry with RPV 
Split Supersymmetry 



would have little impact. In fact, in giving this talk 
in several venues I challenged each audience to pre- 
dict how many of the 21 models considered here were 
constrained by the 13° mixing values or equivalently 
how many evaded constraint. Before answering this 
question, we consider a specific NP example in some 
detail. 

Suppose a vector-like quark of charge Q = -1-2/3 [Tsj 
is added to the SM. Recall that a vector-like quark 
is one whose electric charge is either Q = +2/3 or 
Q — —1/3 and which is an S'[/(2)l singlet. Both 
choices of charge are actually well motivated, as such 
fermions appear explicitly in several NP models. For 
example, weak isosinglets with Q = —1/3 appear in 
Ee GUTs [H, with one for each of the three 
generations {D, S, and B). Weak isosinglets with 
Q = +2/3 occur in Little Higgs theories [H, [l^ in 
which the Standard Model Higgs boson is a pseudo- 
Goldstone boson, and the heavy iso-singlet T quark 
cancels the quadratic divergences generated by the 
top quark in the mass of the Higgs boson. We re- 
strict our attention here to the Q — +2/3 case. Since 
the electroweak quantum number assignments are dif- 
ferent than those for the SM fermions, flavor changing 
neutral current interactions will be generated in the 
left-handed up-quark sector. Thus, there will also be 
FCNC couplings with the Z" boson fl^. These cou- 
plings contain a mixing parameter A„c which is con- 
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FIG. 2: Tree-level contribution from Z^-exchange. 



strained by the unitarity condition 



(15) 



A tree-level contribution to AA/u is thus generated 
from Z^-exchange (see Fig. [2|). It is straightforward 
to calculate that 



.(2/3) 



ri{m,,Mz){D°\ Qi 



BD{Xuc)ri{m,,Mz) (16) 



where we have made use of Eq. Q and Eq. ([7]) . The 
result is displayed in the graph of Fig. [31 which con- 
trasts a ±lcr window (dashed lines) about the HFAG 

central value with the NP prediction x\~, (solid line) , 
as a function of the mixing parameter A„c. The bound 
on A„c from mixing turns out to be roughly two 
orders of magnitude better than that from the CKM 
unitarity constraint. 

Upon performing analogous analyses for the other 
NP models, we arrive at a set of constraints (mainly 
in the form of graphs) like the one depicted in Fig. [31 
It is not possible here to summarize the results for 
all 21 models of Table HVl One must refer to Ref. [2] 
for that. However, we can answer the question raised 
earlier about how many models avoid being meaning- 
fully constrained by the D° mixing data. The answer 




0.01 0.015 0.02 0.025 0.03 0.035 

^uc-lO"" 



FIG. 3: Value of xu as a function of the mixing parameter 
Xuc in units of 10~^ in the Q = -f 2/3 quark singlet model. 
The la experimental bounds are as indicated, with the 
yellow shaded area depicting the region that is excluded. 



is just 4 of the 21 models, which came as a surprise to 
many. These 4 models are Split SUSY, Universal Ex- 
tra Dimensions, LR Symmetric and Flavor-conserving 
Higgs Doublet. 

It is of interest to briefly consider a few of these, in 
order to understand how the D° mixing constraints 
can be evaded: 

1. Split SUSY [131: This is a relatively new vari- 
ant of SUSY (2003-4) in which SUSY breaks at 
a very large scale Ms ^ 1000 TeV. All scalars 
except the Higgs have mass M ^ Ms whereas 
fermions have the usual weak-scale mass. It 
is known that large mixing in SUSY will 
involve squark amplitudes. But since squark 
masses in Split SUSY are huge, the mixing be- 
comes suppressed. 

2. Universal Extra Dimensions [2l|: UED is a vari- 
ant of the idea that TeV~'^-sized extra dimen- 
sions exist. There are no branes appearing in 
this approach, so all SM fields reside in the bulk 
and just one extra dimension is usually consid- 
ered. Each SM field will have an infinity of KK 
excitations. It turns out that GIM cancellations 
suppress all but a few b-quark KK terms, but 
these are CKM inhibited. 

So we see that supressions can arise from more than 
one source and that the suppressing mechanism will 
depend on the specific model. 



V. CONCLUSIONS 

At long last, signals for xd and j/d have been ob- 
served. These experimental findings, although greatly 
welcome, whet our appetite for ever more precise de- 
terminations. Hopefully these will be forthcoming, so 
we can put aside any lingering concerns that all the 
excitement has been the result of statistical fluctua- 
tions. 

The SM analysis, as is so often the case, is not with- 
out its difficulties. At the quark level, theoretical anal- 
ysis in the dimension six sector through NLO gives 
a^D ^ 2/D — 10^^. These values are tiny compared to 
the reported experimental signals. It is evident that 
the triple sum over the operator dimension d„, the 
QCD coupling as and the mass expansion parame- 
ter z of Eq. ^ is slowly convergent. This approach 
remains inconclusive at best. 

A more promising avenue is to study ?/d with the 
hadronic degree of freedom. This yields a plausi- 
ble, and quite possibly correct, explanation for reach- 
ing the ud ~ 0.01 level. Again, however, the ef- 
fect of strong interaction uncertainties mars predictive 
power. 

Finally, the work of Refs. T, '5| has explored which 
NP models can yield sizable values for xd, Ud and 
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which cannot. Charm mixing data has been found ence Foundation under Grant PHY~0555304. 
to infer useful constraints on NP parameters spaces, 
and as should be clear to all, provides a most welcome 
addition to the High Energy Physics community. 
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